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, We consider the density of two-dimensional critical percolation clusters, constrained to touch one 

' or both boundaries, in infinite strips, half-infinite strips, and squares, as well as several related 

quantities for the infinite strip. Our theoretical results follow from conformal field theory, and 
j^', are compared with high-precision numerical simulation. For example, we show that the density 

of clusters touching both boundaries of an infinite strip of unit width (i.e. crossing clusters) is 
proportional to (sin 7ri/)~*/''*{[cos(7r2//2)]^/^ -|- [sin{iiy /2)Y^^ — 1}. We also determine numerically 
, contours for the density of clusters crossing squares and long rectangles with open boundaries on 

' the sides, and compare with theory for the density along an edge. 
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I. INTRODUCTION 



Percolation in two-dimensional systems is an area that has a long history, but remains under very active current 
study. A number of very different methods has been applied to critical 2-D percolation, including conformal field 
I theory (CFT) [1], other field-theoretic methods [2], modular forms [3], computer simulation [4], Stochastic Lowner 
^ . Evolution (SLE) processes [5] and other rigorous methods [6]. (Because the literature is so very extensive we have 
Q ' cited only a few representative works.) 

O , More specifically, there is a great deal of work of recent work studying universal properties of crossing problems 
in critical percolation in two dimensions (i.e., [1, 3, 5, 7, 8, 9, 10, 11, 12]). Another interesting and also practically 
, important universal feature of percolation at the critical point is the density, defined as the number of times a 
J> ■ site belongs to clusters satisfying some specified boundary condition (such as clusters touching certain parts of the 
T-H ] boundary) divided by the total number of samples N , in the limit that N goes to infinity. This problem has been 
. addressed for clusters touching any part of the boundary of a system in various geometries, including rectangles, strips, 
and disks, via conformal field theory [13] and by solving the problem for a Gaussian free field and then transforming 
, to other statistical mechanical models, including percolation [14]. In a recent Letter [30], we considered the problem 
• of clusters simultaneously touching one or two intervals on the boundary of a system, and considered cases where 
\ those intervals shrink to points (anchor points). These results exhibit interesting factorization, arc related to two- 
■ dimensional electrostatics, and highlight the universality of percolation densities. Note that the density at a point 
z of clusters which touch specified parts of the boundary is proportional to the probability of finding a cluster that 
connects those parts of the boundary with a small region around the point z. 

In this paper we consider the problem of the density of critical percolation clusters in various geometries where 
the clusters simultaneously touch both of the boundaries (i.e. crossing clusters), and several related quantities. 
' The first case we consider is an infinite strip, with boundaries parallel to the x-axis at y = and y = 1, so the 
crossing is in the vertical direction. (All our models are defined so that the crossing is vertical. Fig. 4 below illustrates 
the geometries that we consider.) For the infinite strip we find (leaving out an arbitrary normalization constant here 
and elsewhere) that 



Pb(y)==(sin^2/)-^/48 



/ Try \ 1/3 / . Try \ 1/3 
(^eos-j +(sm-j 



(1) 
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This may be compared to a previous result [13, 14] for clusters touching either the upper or lower boundary (or both) 
which is simply given by 

Pe{y) = (sin7ry)-5/48 . (2) 
We also show that the density of clusters touching one boundary irrespective of touching the other is given by 

poiy) = (sin^y)-^/4«(cos^)'^' (3) 

p,{y) ^ (sin^2/)-^/4«(sin^)'^' , (4) 

where po corresponds to those clusters touching the lower boundary at y = 0, and pi corresponds to those clusters 
touching the upper boundary at y = 1. (Note that po is the analog of the order parameter profile {<y)+.f for the Ising 
case (see Eq. (16) in [19].) We also find expressions for clusters touching one boundary and not the other, which are 
combinations of the above results. 

Perhaps the main new theoretical result in the above is (3) (or equivalcntly (4)), which follows straightforwardly 
from the results in [30]. The derivation is given in section II. 

A second type of theoretical prediction gives the density variation along a boundary (the general expressions is in 
(15)). This is used to predict the density along the edge in several geometries (see (17), (18), and (19) below). 

The above theoretical results are found to be consistent with numerical simulations to a high degree of accuracy. 
We include the results of numerical simulations of the density contours of clusters crossing square and rectangular 
systems vertically, with open boundaries on the sides, and compare with theory along the boundaries. 

Our theoretical treatment is related to previous use of conformal field theory to study order parameter profiles in 
various 2-D critical models with edges and similar research [13, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28]. 
(Note that the density which we consider is the expectation value of the spin operator in percolation, which is the 
order parameter in this setting.) Many of these prior results make use of the original Fisher-de Gennes proposal [29] 
for the behavior of the order parameter at criticality near a straight edge. In this paper, we limit ourselves to critical 
percolation. We also include the results of high-precision computer simulations. In addition, the formula (15) for the 
density along the edge of a system is new, to our knowledge. 

Note that, because of the fractal nature of critical percolation clusters, the density of clusters is, strictly speaking, 
zero everywhere in the system. However, if we properly renormalize the density as the lattice mesh size goes to zero, 
the density can remain finite. At the boundaries, for some quantities, this results in the density diverging, as for 
Pe{y) at ?/ = and y = L (but remaining integrable). For pb{y) of equation (1), on the other hand, the renormalized 
density remains finite everywhere. When comparing to numerical simulations, one has to normalize the data so 
that the densities coincide with the theoretical prediction using whatever normalization convention is chosen for the 
theoretical results. The resulting normalization constant, which must be applied to the numerical data, is specific for 
each system and is non-universal. 

In the following sections, we first give the theoretical derivation of our infinite strip formulas above. Then we 
present the numerical results. This is followed by numerical results on square and (long) rectangular systems. These 
are compared with theory for the density along the edges of these systems. We end with a few concluding remarks. 



II. THEORY FOR THE INFINITE STRIP 



We first consider the density of critical percolation clusters which span the sides of an infinite 2-D strip. We can 
find the density predicted by conformal field theory [31] using the results of [30]. In that article we showed that in 
the upper half plane the density p of clusters connected to an interval (.Tq, xj,) is 



p{z,Xa,Xb) CC {Z ~ Z) F[— 

\{Z- Xaj{Xb - Z) J 



(5) 



where the function F[ri) was determined by conformal field theory and takes on one of two forms, 



If we parameterize 2: as z = re'^ and let a;a — > and — > cx), then ?/ = 1 — e^'^ and using (6) we can rewrite (5) as 

p+{r,e,Xa^Q,Xb^oo) cx (rsin6l)-5/''8[cos(6'/2)]i/3 (7) 
P-ir.e^Xa-^Q^Xb^ 00) cx (rsin6l)-5/''S[sin(0/2)]i/3 _ (g) 
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For the positive real axis — > and p+ ^ 0-^/^^ while for the negative real axis — > tt and p_|_ ~ (tt — OY^^'^^. The 
powers here arise from the fixed and free boundary exponents, respectively, in the bulk-boundary operator product 
expansion of the magnetization operator -0 [30]. (More specifically, as it approaches the boundary, -0 ~ 1 or -0 0i,3, 
which have conformal dimensions and 1/3, respectively.) This shows that p+ is the density of clusters attached to 
the positive axis. Because p_ (r, 0) = p+{r, tt — 9), it follows that p- is the density of clusters attached to the negative 
real axis. 

The final density that we need is that of clusters connected arbitrarily to the axis. This is given by {ij{z, z))fixcd oc 
{z - z)-^/-^^ [13, 14, 30] which may also be written 

oc(r sin 0) -5/48 . (9) 

These densities are unnormalized. However for points that arc short distances above the positive (negative) axis, 
but very far from the origin, the relation p+{-) ~ Pa holds. This condition holds since the points arc far from the free 
boundary, and thus dominated by the fixed boundary. It is satisfied by our expressions for and pa, so they 

are properly normalized relative to one another. 

We next map these densities into the infinite strip w Gjx + iy\ x ^ (— oo, oo), y G (0, 1)} using 

w{z) ^ -log{z) . (10) 

TT 

This leads to the expressions for po{y); pi{y), and Pe{y) given by equations (3), (4), and (2), respectively. 
Using these functions we can also determine the densities of clusters that touch one side but not the other, 

Pol(y) - Peiy) - Pi{y) = (sinTry)-^/^^ (l - [cos(7ry/2)]i/3) (H) 

Pio(y) = Peiy) - Poiy) = {smny)-^/'' (l - [sin(^y/2)]i/3) . (12) 

In a similar manner we can find the density of clusters touching both sides, Pb{y)- Adding pQ and pi includes all 
configurations that touch cither side, but double counts the clusters that touch both sides. Subtracting pe leaves only 
those clusters that touch both sides of the strip. Thus 

Pb{y) ^ Po{y) + pi{y) - Peiy) , (13) 

which gives equation (1). 



III. SIMULATIONS FOR THE INFINITE STRIP 



To approximate the infinite strip, we considered rectangular systems with periodic boundary conditions in the 
horizontal direction, for both site and bond percolation on square lattices. Here we report our results for site 
percolation on the square lattice, for a system of 511 (vertical) x 2048 (horizontal) sites, at p = Pc = 0.5927462 
[32]. We generated 500,000 samples to compute the average densities, using a Leath-type of algorithm to find all 
clusters touching the upper and lower boundaries. 

The aspect ratio of the rectangle used in our simulations was 2048/511 = 4.008 . . ., which might seem a bit small. 
However, since the correlation length of the system is governed by the width of the rectangle (511), the effect of the 
finite ratio drops off exponentially with the length of the rectangle, so our results should be very close to those of 
an infinite strip, an expectation which is born out by the results described below. (Furthermore, the probability of 
finding a horizontally wrapping cluster drops exponentially with the aspect ratio, so if a longer rectangle were used, 
very few wrapping clusters would be found and the data would have poor statistics.) 

The agreement of predicted and simulated values is excellent. In Fig. 1 we show the results for pbiy) with no 
adjustments, other than an overall normalization (in particular, the extrapolation length a, described below, is set to 
zero). Fig. 2 (upper set of points) shows the ratio of simulation to theoretical results, equation (13). Most points fall 
within 1%, except those right near the boundary. However, we can do better. 

When comparing simulations on a (necessarily) finite lattice with the results of a continuum theory, there is the 
question, due to finite-size and edge effects, of what value of the continuous variable y should correspond to the lattice 
variable Y — specifically, where the boundary of the system should be placed. On the lattice, the density goes to zero 
at y = and Y ~ 512. A naive assignment of the continuum coordinate would therefore he y ~ Y/512. However, 
we can get a better fit to the data near the boundaries by assuming that the effective boundary is a distance a (in 
units of the lattice spacing) beyond the lattice boundaries — that is, at y = —a and Y = 512 + a. The distance a is 
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FIG. 1: (Color online). 

Pb{y) vs. y. Open violet circles, theoretical values from equation (13) (here normalized to 1 at the center point), using (14) 

with a = 0. Blue dots: results from simulations. 



an effective "extrapolation length" where the continuum density far from the wall extrapolates to the value zero [33] . 
This is accomplished by defining the continuum variable y by 

y = (y + a)/(512 + 2a) . (14) 

Then, Y = —a corresponds to y = and Y = 512 + a corresponds to y = 1. Note, Y ^ corresponds to 
y = a/(512 + 2a) and Y = 512 corresponds to y = (512 + a)/(512 + 2a) = 1 - a/(512 + 2a), and so the theoretical 
extrapolated density pb{y) will be greater than zero at these points on the actual boundary. The spacing between all 
points is stretched by a small amount because of the denominator in equation (14), but this stretching does not have 
much effect on the behavior of ph near the center. The main effect of a on the shape of the theoretical curves of pb is 
near the boundaries. 

By choosing an extrapolation length of a = 0.26, we can get a much better fit of the data, as can be seen in Fig. 2 
(lower set of points) . A plot of the data analogous to Fig. 1 puts most of the data points right in the center of the 
theoretical circles, but would barely be visible when plotted on this scale. With a = 0.26, the error is now reduced to 
less than 0.1%, except right near the boundaries, as can be seen in Fig. 2. A more thorough study of the extrapolation 
coefficient a would require the study of different sized lattices, and the demonstration that a is independent of the 
lattice size. We have not carried this out. Note, however, that a constant a implies that if one keeps the physical size 
of the lattice fixed (so that the increasing number of lattice points makes the mesh size go to zero), the extrapolation 
length, measured in physical units, will also go to zero. 

Note also that the distance in the y directions was 511 rather than 512 because we used one row at y = 0, in 
conjunction with the periodic boundary conditions, to represent both horizontal boundaries of the system on the 
lattice. That is, F = 1 and Y = 511 are the lowest and highest rows where we occupy sites in the system, where Y 
represents the lattice coordinate in the vertical direction. 

We note that these simulations were carried out before the theoretical predictions were made. 

We have also carried out simulations measuring the density of clusters touching one edge, po{y)- The results are 
shown in Fig. 3. We also plot the results of the theoretical prediction, equation (3), on the same plot, and find 
agreement within 0.5% without using an extrapolation length a. 
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FIG. 2: (Color online). Ratio of simulation to theoretical results for Pb{y) with a — (upper set of points) and a — 0.26 (lower 
set of points) 




FIG. 3: (Color online). Density of clusters touching lower boundary, po{y), as a function of y, both simulation (dots) and 
theory (open circles), equation (3). 
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FIG. 4; Sketches of the cases considered. Solid (dashed) boundary lines represent fixed (open) boundary conditions. Curved 
lines indicate crossing clusters; the density pb is evaluated along the lines with arrowheads, (a) is the infinite strip, cf. e.g. 
equation (f); (b) the square (Fig. 6), (c) and (d) vertical half-infinite strips ((17) and (18), respectively); (e) horizontal 
half- infinite strip (19). 



IV. PERCOLATION ON A SQUARE AND SEMI-INFINITE STRIP 

In this section, we consider the density of crossing clusters on a square with open boundaries and also on a (long) 
rectangle. We compare the numerical results with the predictions of conformal field theory for the density along an 
edge. The various cases considered are illustrated in Fig. 4. 

Note that, as mentioned, the crossing is always in the vertical direction. In a slight abuse of notation, we use pb for 
the density of a clusters that touch both anchoring intervals in all cases. The different situations may be distinguished 
by the arguments of pb, e.g. pb{y) for the infinite strip, where there is no x dependence, and pb{x,y — 0) along the 
bottom of the semi-infinite strip as in (17) below. 

Our simulations of percolation densities on an open square examined clusters that cross in the vertical direction. 
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FIG. 5: (Color online) Contours of constant densities 0.625, 0.75, 0.875, 15/16=0.9375, 31/32=0.96875, 63/64 = 0.984375, 
127/128 = 0.9921875, and 1 (outside to center) of clusters touching both the top and bottom edges, with open b.c. on the sides, 
for a system of 511 x 511 sites. 



with open boundary conditions on the left- and right-hand sides. We considered site percolation on a square lattice 
of 511 X 511, with 2,000,000 samples generated. The resulting contours are shown in Fig. 5. As in the periodic case, 
the density goes to zero at the upper and lower boundaries because, compared to an infinite system, these boundaries 
intersect many possible crossing paths, leading to large holes in the clusters near the boundaries. As a consequence, 
relatively few points on the boundary will be part of the crossing clusters, and in the limit that the mesh size goes to 
zero, their density evidently goes to zero. The density also goes to zero on the sides because of the open conditions 
there. Interestingly, the contour curves are almost symmetric in the horizontal and vertical directions, indicating that 
the anchoring and open boundaries have a similar effect on the density. 

We have not carried out a field-theory calculation to find the density inside the square or rectangular systems. To 
do so requires a six-point function whose calculation would be unwieldy. 

It is however relatively easy to calculate the variation of the density along the bottom edge of the square, now 
normalized so that the density remains non-zero as the mesh size goes to zero. To do this, we consider the density of 
crossings from one of the anchoring intervals to a single point on the other interval, using the boundary spin operator. 
Now crossing from the top edge to one point x on the bottom edge (which is given by a three-point function, depending 
only on x and A) automatically implies crossing from the top to bottom (which is given by a four-point function, 
depending only on A). Therefore the density at x is proportional to the ratio of the former to the latter. It follows 
generally that, up to a (A-dependent) multiplicative constant, one has [34] 



Here z(w) maps the region of interest {w) onto the 1/2-plane (z), x is the i/j-coordinate along the anchoring interval 
of interest, and A is the conformally invariant cross-ratio for the anchoring points. For a rectangle, this depends on 
the aspect ratio r = A'(\/l - A)/A'(VA) [1, 35]. 
The mapping for the square is 




(15) 




(16) 
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with p(u; (72,53) the Weicrstrass elhptic function and K{z) the eUiptic integral function. This mapping takes a unit 
square x,y € (0,1) into the upper half plane. For the square A = 1/2, and we can take x G (0,1). In Fig. 6 we 
compare the measurement and theory. Clearly the agreement is excellent. 

In the case of a half-infinite strip x G (0, l),y € (0, 00), the density of sites along the a;-axis belonging to clusters 
crossing vertically is found from (15) using z^w) = sm^{nw/2), and A = 0. This gives 

Pb(x,y = 0) = (sin^x)i/3 . (17) 

Of course, for an infinite strip, the probability of crossing (in the long direction) is zero, so one must consider the 
limit of a large system, and calculate the density given that crossing takes place, and take the limit that the length of 
the rectangle goes to infinity. It turns out that numerically, the density at the edge for the square, equations (15,16), 
differs only very slightly from that of the half-infinite strip, given by equation (17). From the point of view of the 
density along an anchoring edge, the square is not much different from the half-infinite strip. 

For J/ ^ in the above half-infinite strip (or equivalently for any y for a fully infinite strip in the vertical direction) 
one can also find the density along the x-direction of the vertically crossing clusters in closed form 

Ph(x,y>0) = (sinTra;)"/^^ . (18) 

This function may be found by transforming the density of clusters connecting two boundary points, derived in [30], 
into the infinite strip. We then take the limit as the two anchoring points move infinitely far away in opposite 
directions, while normalizing the density so that it remains finite. (Related order-parameter profiles for the Ising 
case are studied in [22].) Interestingly, a plot of this density profile (written as a function of y rather than x) is 
very similar in appearance to that of the vertically crossing clusters pb{y) given in equation (1) and plotted in Fig. 1. 
When normalized so that they have the same value at y = 1/2, the maximum difference between the two is at y = 
and y = 1, where (18) is only 1.5 % below (1). This small difference indicates that open boundaries and anchoring 
boundaries have similar effects on the density of the crossing percolation clusters, and is consistent with the near 
symmetry seen in the contours in Fig. 5. 

We can also find the density along the left and right (open) edges. For the case of a half-infinite strip rotated by 
90° with respect to the one above, (y € (0, l),a; € (0,oo)), we find 




which is similar in form to pb{y) of equation (1) (which in fact corresponds to x ^ for the geometry considered here) 
but with different exponents. This similarity arises because the derivations of (19) and (13) are virtually identical, 
except that for (19) we leave a free interval between the two anchoring intervals (where the boundary spin operator 
sits) which is mapped to the end of the half-infinite strip using sine functions. Comparison with numerical data 
(not shown) for the density at the short edge of a rectangle of aspect ratio 8 to approximate the infinite strip shows 
excellent agreement with (19). 



V. FURTHER COMMENTS AND CONCLUSIONS 



If we consider densities raised to the sixth power (compare Rcf. [30]), we find a Pythagorean-like relation involving 
Pe{y)^ 1 Po(]jY , and pi(y)^ (which we present without interpretation), 

P^{yf + Pi{yf ^ Pe{yf ■ (20) 

There seems to be no simple relation involving pb{y) other than its basic definition given in equation (13). For the 
corresponding quantities at the edge of the strip as in (19), the power in equation (20) is 3 rather than 6. 

Although the overall normalization of a density such as Pb{y) (for the infinite strip, see (1)) is arbitrary, we can fix 
its value by requiring that 



Pb{y)dy = 1 



(21) 



That is, we define pb{y) = B sin(7ry) ^/^®(cos(7ry/2)^/^ + sin(7ry/2)^/'^ — 1) where i? is a constant. Then equation (21) 
yields 



B = 



, 32r(i)r(i) 



V^r(i) 0Fr(3L)r(i) 

1.46408902 .... 



(22) 
(23) 
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FIG. 6: (Color online) Density of vertically crossing clusters pb{x,0) along the lower boundary, in a square system with open 
boundaries on the horizontal sides. Red line: theory, equation (15), black dots: simulation results. 



Another choice of B is to make Pb(l/2) = 1, which yields B = (2^/^ - l)-i = 1.27910371 . . . 

In many problems of percolation density profiles, the density goes to infinity at a boundary point, such as occurs for 
Po{y) when y 0. Interestingly, in the case considered here of clusters touching both boundaries, Pb{y), the density 
goes to zero at those boundaries and remains finite everywhere. 

To highlight the difference between densities of all clusters touching a boundary vs. the densities of crossing clusters 
touching one boundary, we consider the limit that the strip width becomes infinite, so that the system becomes a 
half-plane. Because we have written our results for a strip of fixed (unit) width, this density is given by the behavior 
of p for small y. The density of all clusters touching the x-axis is thus found from equation (3) (or [13, 14]) to be 

Po - y-'"^ , (24) 

where we have left off the coefficient because the normalization is arbitrary. This density diverges at y = because 
it is much more likely to find sites belonging to these clusters near the .T-axis. On the other hand, the density of 
crossing clusters touching the x-axis is found from equation (1) in the limit ?/ — > 0, 

Pb - y''"' ■ (25) 

In this case, the density increases as y increases, in contrast to (24), and goes to zero at the anchoring boundary, for 
the reason mentioned above. 

Behavior identical to (25) also follows from the small- a; expansion of pb{x,y ^ 0) given by (18), which represents 
the behavior of the density of the vertically spanning clusters at the open boundary. Thus, near the boundaries (but 
away from the corners), the open and anchoring boundary conditions have identical effects on the density of the 
crossing clusters. 

In conclusion, we have studied the density of vertically percolating clusters in a square system, as well as for half- 
infinite and infinite strips extending in either the horizontal or vertical direction. The various cases considered are 
illustrated in Fig. 4. 

For the infinite strip. Fig. 4(a), the density for crossing clusters is given by (19) and compared with numerical 
simulations in Figs. 1 and 2. For the square, Fig. 4(b), we find theoretical results for the anchoring edge densities 
(see (15) and (16) and Fig. 6), as well as numerical results for the density in the interior (Fig. 5), which exhibit an 
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interesting near symmetry. For the half-infinite strip in the horizontal direction, Fig. 4(e), the density pb{x = 0, y) at 
the left open boundary is given by (19), and at a; ^ (or equivalently, for an infinite strip), the density is given by 
(1). For a half-infinite strip in the vertical direction, the density along the lower anchoring boundary (Fig. 4(c)) is 
given by (17) while for y 3> (Fig. 4(d)) the density is given by (18). For the half-infinite systems, the densities near 
a wall are given by the same power-law (25) regardless of whether it anchors the crossing clusters or is open. Note 
that all of our theoretical predictions were confirmed by computer simulation. 

For the future, it would be interesting to study analogous properties for Fortuin-Kastelcyn (FK) clusters of the 
critical Ising and Potts models. 
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